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SUMMARY 

A graphical method is presented for calculating the linearized 
lift distribution on thin wings at supersonic speeds. The technique 
may be applied to all wing regions except those influenced by inter- 
acting flow fields off the wing-plan boundaries. The lifting- 
pressure coefficients are obtained as the sum of a graphical line 
integration and several terms that are functions of only the plan 
form. 


Pressure coefficients obtained by the graphical method for a 
swept trapezoidal wing and for a wing with straight- swept leading 
edge and parabolic wing-tip are compared with those derived by 
closed-form integration. The lift distributions for two plan forms 
previously' unsolved are also included. 


INTRODUCTION 

The evaluation of the theoretical performance of arbitrarily 
shaped wings at supersonic speeds has been retarded by the complex- 
ities involved in analysis when the wing boundaries are other than 
certain prescribed curves, usually straight lines. Solutions for 
the velocity potential on an arbitrary wing surface are presented 
in references 1 and 2 except for those regions influenced by Inter- 
acting flow fields off the wing-plan boundaries; the integrals 
involved, however, have been tabulated only for special plan forms . 
Although the integrals may be so reduced that numerical evaluation 
of the velocity potential is feasible, the procedure for obtaining 
the lift distribution requires an additional numerical partial dif- 
ferentiation that is cumbersome and of doubtful accuracy. 

The lift distributions for a family of thin wings may be obtained 
by the methods of reference 2 as the sum of an algebraic function 
and a line integral. The line integral may be evaluated by numerical 
or graphical methods. 
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A gr aph! cal method, formulated at the NACA Cleveland laboratory 
for evaluating the lift distribution for wings of arbitrary plan 
form within the limitations of reference 2 is presented. Also 
included for a Mach number of V1T are curves showing the relations 
between the algebraic functions and the geometry of the wing. The 
application of the method to the various types of region encountered 
in analyses of general wings is discussed. As illustrations, the 
pressure coefficients obtained by this graphical method are compared 
with those derived by closed-form integration for a swept trapezoidal 
wing for a wing with a straight leading edge and parabolic wing- 
tip boundaries. Lift distributions of two other plan forms are also 
included to illustrate the application of the method. 


SYMBOLS 

The following symbols and subscripts are used in this report: 


& 2 , ag, , . coefficients in series expansion defining curva- 

ture of leading edge 


a*, b o', d f quantities used in determining factor AFg 


C 1 • • • C n coefficients in series representation of AFg 

Cp pressure coefficient 


AC t 


contribution to Cp of portion of leading edge 
included in segment of width € 


*1 


AF 2 


factor representing integral-function contribu- 
tion to Cp 

factor representing contribution of integration 
over width e to Cp 

error in Fg due ‘ to neglecting ourvature 


H 


quantity representing numerical integration 
along leading edge, ArjSb 


h 


parameter defining equilateral hyperbola 


k l 


constant determining sweep cf straight leading 
edge 
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M 

m 

Q(0 

B 

r 

U 

u, v 

UgCv) or v 2 (u) 
v x (u) or u x (v) 

x, y 
y* 

a 

P 

e 

6 1 > e 2 
q> 

Subscripts: 

0 

L 


free- stream Macb number 
slope of wing edge, 

series defining curvature of leading edge 

v w - Tl (u 2 ) 
ratio of distances, — 

local radius of curvature 

free-stream velocity (parallel to x-axis) 

oblique coordinates whose axes lie parallel to 
Lfoch lines 

equation of tip and trailing edge of wing 

equation of supersonic leading edge of wing 
(inboard of foremost tangent Mach wave) 

Cartesian coordinates 

transformed y- coordinate 

angle of attack, radians 

Vm 2 - 1 

width of strip chosen in evaluating F 2 

local angles between wing edge and free-stream 
flow direction 

Cartesian coor dina tes of point sources 
perturbation-velocity potential 

intersection of grid edge with wing-plan boundary 
leading edge of wing 
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l left 

r right 

w wing 


ANALYSIS 

As a result of the linearized theory, the effeot of angle of 
attack on the lift distribution of thin wings at supersonic speeds 
depends upon the plan-form boundaries and is independent of the 
wing-section slopes. Determination of the pressure distribution of 
a thin flat plate therefore solves the problem of the lift distri- 
bution of the finite- thickness wing. Reference 2 shows that the 
local pressure coefficients of several classes of wings with arbi- 
trary plan boundaries may be evaluated by line integrals of explicit 
functions. 

A graphical method of evaluating the lift-distribution functions 
presented in reference 2 is developed herein. In order to illus- 
trate the method developed, a simple wing including regions influ- 
enced by a supersonic leading edge, a subsonic leading edge, and a 
subsonic trailing edge (fig. l) is discussed in detail. The essen- 
tial elements in the calculations for each region are included in 
the analysis for the region B influenced by the subsonic leading 
edge. 

The pressure coefficient for region B as derived in refer- 
ence 2 is 

o r 

+ 1 i 1 - 

where U2 is evaluated at v = v w . Except for the symbol tj, 
equation (l) is expressed in a set of oblique coordinates (fig. 2) 
whose axes are parallel to the Mach lines. The transformation 
equations relating the oblique and Cartesian coordinates are 
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co 

ro 


u= 

V 

-§ « + ^ 

| = I (T + tt) 

n 

= i (T - u) 


T w 

- ft + to) 

* - I ( T v + 

y 

= M 




> 


J 


( 2 ) 


The significance of the terms of equation (l) is clarified in 
figure 3. The line integral is to "be evaluated along the portion 
of the supersonic leading edge from a to b. The quantity 
[v w - V!(u 2 )] can be interpreted as distance be and the quantity 
(uy. - u 2 ) can be interpreted as the distance cP. The ratio of 
these distances is called E. For a given Mach number, the quan- 


tity ^ 1 - depends upon only the slope of the wing boundary 

at the point (u 2 ,v w ). 


By use of equations (2), the second term of equation (l) may 
be eapressed in terms of E and the angle that the wing edge 
makes with the free-stream flow direction (parallel to the x-axls) 
as 



dT v/ 


Vi 


ifi 


tan 9i 


M (p tan 8]_ + 1) 


Vi 


( 3 ) 


where 0p is positive counterclockwise . This espression, eval- 
uated for M = 4%, appears in figure 4 . 

The quantities (u^. - u) and (v v - v) in the line integral 
of equation (l) are the coordinate distances (fig. 3) of each 
element dq to the point (u w .,v w ). The curves of 


1 

V(uw - u) (y v 




= • constant 


= h 
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in the u,v coordinate system are a family of equilateral hyper- 
bolas asymptotic to the lines u = v = v w . A family of these 
hyperbolas is shorn in figure 5. 

If a "grid" of this type is placed with its origin at (u^v^) 
and so alined that its axes are parallel to the u and v axes, 

the value of ■■■■ _ ■ _ - ~ . for any element dq (or dy) is 

V(«w ~ u ) ( v w “ T ) 

the constant h for the hyperbola passing through that element. 

The Siam of the values of h taken near the center of each element 
for equal increments drj along the wing boundary then evaluates 
the line integral of equation (l). The superposition of a wing- 
plan 'boundary on a hyperbolic grid for M = Vt' is shown in fig- 
ure 6. A constant value of Ar| is marked on the grid for conven- 
ience in summation. A part of the grid employed for calculations 
at M = V2 appears in figure 7 . 

As shown in figure 5, the value of h increases to infinity 
as the edge of the grid is approached. In order to perform a 
numerical integration, the infinities must be excluded. Therefore 
the numerical summation should not start with element 1 (fig. 6) 
but should skip a distance large enough to enable the element 
value to be easily determined. The part of the Integral omitted 
must be compensated by an analytical expression that represents 
the contribution of a section of the leading edge aa’ included 
in a strip of width e (fig. 6). 

The contribution to the line integral of the strip of width e 
may be evaluated by assuming the wing boundary to be a straight 
line near u = u„.. Prom the derivation presented in appendix A, 
this contribution is 


1 
2 

( 4 ) 

dy 

where m = the tangent of the angle (positive counter- 

clockwise) that the leading edge makes, with the flow direction at 
the intersection with the edge of the grid (point a, fig. 6), and 
(v w - Yq) is the distance aP from that point to the grid origin 

(fig. 6). The value of Fg for various values of 62 , (v v - v 0 ). 




M V( m P - + 1) 


tan' 


-1 


2f3m(v w - Vq) 
M e (1 + pm) 


- 1 
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and € is shown in figure 8 for M « V2. The dashed line in fig- 
ures 8(h) and 8(c) represents the limit for which the arc-tangent 
term in equation (4) is real. 

Jor the subsonic leading edge (region B, fig. l), equations (l), 
(3), and (4) may he combined and the expression for the pressure 
coefficient becomes 


Op - - § ( H + Tl + J? 2 ) (S) 

where the numerical summation along the leading edge is repre- 
sented by 

H = Aq (h-L + h 2 + h 3 + . . .) = AqSh 

If nonuniform increments of q are desired, a weighting of each 
value of h is required. 

The analysis of the region of a wing tip influenced by a sub- 
sonic trailing edge (region C, fig. l) may be appreciably shortened 
if the Kutta-Joukowski condition Is imposed. In reference 2, it 
is shown that the solution which satisfies the Kutta-Joukowski con- 
dition does not contain the term Fq. Therefore Fq Is taken as 
zero In this region and the remaining evaluation is unchanged. 

When evaluating C p in a region where v<0 (region A, 
fig. 1, or fig. 9), the value Fq does not exist because no sub- 
sonic leading edge is included in the forward Mach cone from 
point P. The integration must be conducted over the entire wing 
leading edge between the limits of the grid a and b. Because 
h becomes infinite at the right as well as at the left limit of 
the grid, a part of the integration at the right must be replaced 
by an analytical expression similar to equation (4). Equation (4) 
can be shown to hold at this limit if (uy- - u q ) is substituted 
for (v w - vq) and ( 180 ° - 6q) is substituted for 0 2 . If the 
value of this expression is designated F 2 r the pressure coef- 
ficient in the region where v<0 becomes* 

C p = - ^ (H + F 2 + F 2 , r ) (5a) 

In general, a complete wing in supersonic flow may have regions 
under the influence of both wing tips. A summary of the method of 
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obtaining lift distributions in the various types of flow field 
commonly encountered is given in appendix B. Also included is a 
numerical evaluation of the pressure coefficient at point P of 
figure 10. 

The effect of Mach number on lift distribution for a given 
plan form may be determined by two alternate procedures. One 
method is to construct an integration grid and curves of and 

F 2 for the Mach number desired. If a large number of analyses at 
a given Mach number are to be made, this method may be preferable. 
If only a very few analyses at a given Mach number are to be made, 
however, the grid and curves of Fj_ and F 2 ma y be used as con- 
structed for M = V 2 , provided that compensating corrections are 
made. These corrections consist in: (l) so transforming the wing 

boundary that y* = By; and (2) analyzing the wing as at M = 
and dividing the resultant Cp/a by the value of 3. 

Some considerations' that simplify the application of the 
method appear in appendix C. Construction of the integration grid, 
choice of Atj and e, drawing size, and a special example of the 
subsonic trailing edge are discussed. 


ACCURACY OF METHOD 

In order to determine the accuracy of the method, the pressure- 
coefficient at points on two wings for which analytical expressions 
were obtainable were graphically computed and compared with the 
analytical values. The wings and the points considered appear in 
figure 11. The leading edges of both wings are swept back 30° and 
are the same except that wing A has a straight tip and wing B has 
a parabolic tip passing through similar points on the leading and 
trailing edges. The points considered are at the same x,y values 
for each wing. The results of the computations for Cp are pre- 
sented in table I, 

The magnitude of the personal error incurred in application 
of the method was evaluated by trial. Values of Cp/a were inde- 
pendently obtained by nonprofessional personnel using the graphical 
method, and were also analytically calculated. For wing A, the results 
of one of the computers contained an average error of 0.46 percent 
and a maximum error of 0.71 percent (table I). Another computer 
averaged an error of 0.53 percent with a maximum error of 1.52 per- 
cent. Results of wing B gave about the same accuracy (table I). 

These results were obtained from drawings with a wing chord of 
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4g Inches. The computers estimated a time rate of approximately 

5 minutes per point when a large number of points were being 
analyzed. 

Inasmuch as the first-order expression for F 2 (which also 
includes higher-order terms) assumed a straight leading edge, this 
expression was exact for the wings shown in figure 11. The accuracy 
of this expression for a curved edge has been investigated and is 
shown in appendix A. For a parabolic leading edge, the error AF 2 
in F 2 (including fourth-order terms) as a function of the ratio 
c/ro (where tq is the radius of curvature) is presented in fig- 
ure 12. If e/rQ is small, the error in Fg due to curvature 
will be small except when 0 2 approaches the Mach angle, as indi- 
cated in figure 12. The corresponding error in Cp will, of course, 
be a smaller percentage. A similar trend probably holds for edges 
of different curvature. 


APPLICATIONS 

As examples of the method, the lift distributions 2Cp/a of a 
circular-plan-form wing and a straight wing with a circular tip are 
shown in figures 13 and 14, respectively, at M 

For the circular wing, in the region inboard of station D, the 
pressure at a given chordwise station increases in the outboard 
direction. For a given spanwlse station, the pressure is a maximum 
at the leading edge and drops off in the flow direction. The rate 
of this decrease is a minimum at the wing root. 

Lines of constant pressure for the circular-tip wing swept 
back 30° at a Ifech number of V2 are shown in figure 14. The 
Kutta-Joukowski condition was assumed for analysis of regions 
influenced by the subsonic trailing edge. The only region of high 
pressures and high pressure gradients is a small leading-edge 
region between the tangency of the foremost Mach wave and the start 
of the subsonic trailing edge. The pressure gradient over most of 
the tip is roughly in the spanwise direction. 


CONCLUDING DISCUSSION 

A graphical method based on linearized supersonic-flow theory 
has been developed for calculating the lift distributions on thin 
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wings at supers onio speeds. The method is applicable to all of the 
regions of arbitrary wings except those affected by interacting 
flow fields off the wing-plan boundaries. The determination of 
the lifting-pressure coefficient for a given point has been, found 
to require 5 to 10 minutes with a resultant average error of less 
than 1 percent as determined by comparison with known analytical 
solutions. 


Plight Propulsion Be search Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, May 13, 1948. 
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APPENDIX A 

w CONTRIBUTION TO C p OF LEADING-EDGE SEGMENT 

^ INCLUDED IN STRIP OF WIDTH e 

In order to determine the contribution to C p of a segment of 
width €, the integral of equation (1) must be evaluated. In the 
x,y coordinate system, this integral is 


— / dr) 

31 J V(* - |) 2 -P 2 (y - tj ) 2 


The integral is considered in the region near 

(x - i Q ) 2 -p 2 (y - i) 0 ) 2 = 0 


(Al) 


(A2) 


The subscript 0 refers to the point of intersection of the grid 
edge and the wing-plan boundary. In particular, the positive root 
that defines the left forward Mach line from the point fc,y) is 


(x - | 0 ) - 0 (y - no) = | ( v w - T o) 


(A3) 


Treatment of e as a variable in the vicinity of the point 
0) 0 # I 0 ) yields 


where 


and 


t] = t } 0 + e 

| =| 0 + | + 3(0 




m = 




Q(e) = &2 e2 + a 3 e3 + »4£ 4 + . . . 


J 


(A4) 
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Substitution of equation (A4) in equation (Al) gives 



Now, if Q(e) = 0, tbe leading edge is straight and the inte- 
gral may be written 




This function, which is equation (4) of the text, is shown in fig- 
ure 8 for various values of an(i ( T w " Y o ) / € ■ for M “ ^2* 

If the leading edge is curved, Q(c) ^ 0 and equation (AS) 
must be rearranged. The quantity in the radical may be expressed 
as 
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b* =2 (X - ! 0 ) Q - 2 ('I) Q - Q 2 
2 (x - | 0 )(mp - 1) 



Then equation (AS) may he written 



Expanded in a power series, 



Equation (AS) then he comes 



(ASa) 

The first term is identical with the integral of equation (A6). 

The contribution of the curvature to Fg nay therefore he written 
as 
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(A9) 


*k * 1 

The expressions for —7 and -==r In equation (AS) may "be 

a Va ’ 

expanded In a paver series to yield 

aT = - t 0 ) a 2] c + ^ x - ^o) a 3 “ 5 " fr ^ 0^ a 2 J * Z 

+ [(x - | 0 )a 4 - ^ a 2 2 - (fr) ( x - £ 0 ) a 3 + (§7) “ 



_1 _ 1 _ 
VT V5 7 


1 





where 


o' 


-(I + mg) 
2 m (x - | 0 ) 


(A10) 


Upon substitution of equation (A10) in equation (A9) and multi- 
plication and collection of terms according to powers of e, there 
results 
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A3?«- 


23 (V 
~ M \o\ 


|rv i 2 2 

‘ SC-^ 2 + 5C 2 e 2 + 7C 3 € 2 + • • 



de 


_ ifi /Li\ 2 ( c -2 , n -2 , n e 2 . 
Id [ci 7 ) V ° ie Gz€ ° 3 

■where 

Oi - | U - 1 0 ) a 2 


°2 'S 


°s- 7 l 


1 I (33m - l)a 2 


3ap 2 m 


a _ + 

4 (3m - 1)_ 

(x T Sq) 


(A9a) 


4m 


3(3m + l) (53m - 3)a 2 (73m - l)a 2 (33m - l)a 3 

— + + - 


32m 2 (x - t Q ) 


16(3m - 1) 


4m 


3a 2 a3m 5a 2 m 


2(3m - 1 ) 8 ( pn _ D 2 J 


(x - S 0 ) 


(All) 


C 4 = 


The series of equation (A9a) converges as long as h'/a' and 
d'e/c' (defined "by equations (AS) and (A10)) are less than unity. 
For rapid convergence, these ratios should "both he small. 

As a particular example (taken at M = V ¥), the leading-edge 
curve may he represented hy the parabola 


S « So + £ + a 2 e2 

221 

The coefficients of equation (All) evaluated in terms of the oblique 
coordinate distances then become 
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C4 = • • • 

* 

When the values of the coefficients are substituted In equation (A9a), 
AFg results. 

For a parabola with a given initial slope m at the point 
(Iq, tj 0 ), the quantity ag controls the radius of curvature r 0 
at that point. The ratio AFg/Fg represents the relative error 
resulting from the assumption of a straight leading edge for a 
width e when the actual edge is curved. This value from equa- 
tions (A7) and (A9a) has been plotted against the ratio e/r 0 * in 
figure 12. 
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APPENDIX B 


TEEATMENT OP COMPLETE WINGS 

A tabular summation of the various flow fields encountered in 
wing problems is presented in table II and the graphical solution 
in each region is indicated. In the reference sketches for 
M = V2, P is the point at which the pressure coefficient is to be 
evaluated. The numerical line integration along the supersonic 
leading edge is performed from point a to point b. The direc- 
tion of integration of the line integrals is indicated by arrow- 
heads on heavy lines. For a subsonic trailing edge (sketches C, F, 

G, H), the Kutta-Joukowski condition was imposed and made F^ = 0 
for that edge. The faotor Fg vanishes when the limit of the 
leading-edge line integral is other than the integration-grid edge 
(sketches D, E, F, G } H), because no infinities are then encountered. 
When the forward-reflected Mach lines from the point P cross on 
the wing surface (sketches E, F), the integration E is considered 
negative. 

The occurrence of these flow fields on a schematic wing is 
illustrated in figure IS. Each flow region is identified by the 
corresponding letter from table II. In this emmpls, flow 
region H does not occur. 

As an example of the method, ■ numerical calculations of the 
pressure coefficient at point P for the wing boundary shown in 
figure 10 are 

6 1 = 7.5° 
be = 0.0419 
cd = 0.084 


B = ££ = 0.4988 
cd 




18 


NACA TN No. 1676 


9 2 - 125.5° 

(v w - Vq) = 0.0772 
C = 0.01 

H = 0.01 [(0.4) (31) + 25.0 + 21.0 + 18.85 + 17.55 
+ 16.75 + 16.50 + 16.6] = 1.446 

. (H + F-, + F ? )= 1.048 

art j. c 
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APPENDIX C 


SIMELIPYING- CONSIDERATIONS IN APPLICATION OP METHOD 

A simple method of obtaining coordinates of hyperbolas to con- 
struct an integration grid employs log- log graph paper. If a 
straight line of slope equal to -1 is drawn on this paper, the 
coordinates of this line, when plotted in Cartesian coordinates, 
describe the desired equilateral hyperbola. 

The addition to the grid of lines of constant v and u 
(parallel to the grid edges) with accompanying scales at the grid 
edges increases the rapidity of evaluation, because all signifi- 
cant distances may then be read directly. A part of the grid 
employed is given in figure 7 . 

In choosing the value of e , the first consideration Is that 
€ be large enough that the value of h may be easily determined 
for the first subsequent element. It is also seen from figure 8 
that if e is taken as a simple decimal such as 0.01 or 0.02, and 
(v w - vq) is known, the ratio (v w - Vq) / e may be rapidly computed 
and hence Pg is rapidly determined. Prom figure 6, point a 
will not, in general, fall at an integral multiple of Aq$ then, 
if 6 = NAt) where N is an integer, point a' will not, in 
general, fall on an integral multiple of Aq. Thus a section whose 
width is less than At) appears before the first complete leading- 
edge element and remains to be Included in the numerical Integra- 
tion. If this small element is evaluated in the normal manner and 
weighted according to its width, no discrepancy will occur. This 
weighting process is illustrated by the example presented in 
appendix B. 

The most prevalent personal error in the application of the 
method was in the measurement of angles ©3. and ©g, because 
on small drawings it is difficult to Judge accurately the tangency 
of a curve. Large-scale drawings are therefore recommended for 
measuring the angles. These drawings may also be used for the 
integration, although the increase in acouracy may not Justify the 
accompanying inconvenience. 

Por a straight supersonic leading edge, equation ( 1 ) is 
integrable and the expression for Cp becomes 
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T w + k l u 2 
- u 2 


k l + 1 -l 
tan 




'^lKr - u 2> 

y v + k l u 2 


(Cl) 

In the region affected by the subsonic trailing edge (region abc 
of ‘Figs. 16 and 17) the first term in equation (Cl) vanishes if 
the Kutta-Joukowski condition is imposed. The quantity in the 
radical of th e arc-tangent term is proportional to the ratio of 
distances sP/rs (fig. 16) and hence any wing composed of the 
same leading edge and a straight line through point a would give 
the same value of Cp for any point on line sP. If lines of 

constant pressure are desired, the arc-tangent term in equa- 
tion (Cl) must be differentiated and set equal to zero. The result 
involves a derivative du 2 /dv w , which is determined by the slope 
of the edge at point s. Equivalent lines of constant pressure 
along line sP may thus be obtained by replacing the tip by a 
straight-line tip that is tangent at point s. The pressure along 
any v = constant line is therefore conical about a point deter- 
mined by the intersection of the extension of the leading edge and 
the tangent to the tip at u 2 (v w ). Along v = v x and v = v 2 
(fig. 17), the constant-pressure lines are conical about points t^ 
and t 2 . By this relation, the constant-pressure lines may be 
quickly constructed and evaluation is necessary only along one line 
of constant x to determine the pressure field. The relation holds 
only in the region influenced by the subsonic trailing edge when 
the Kutta-Joukowski condition applies, for differentiation of 
equation (Cl) with the first term included yields a second deriv- 
ative that generally cannot be evaluated from straight-line 
relations . 
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TABLE I 

VALUES OF Cp/a OBTAINED FOR WINGS A AND B 


Station 

u 

V 

Wing A 

Wing B 

C p /a 

Error 
(percent ) 

Cp/a 

Error 

(percent) 

Analytic 

Graphic 

Analytic 

Graphic 

1 

8 

0 

2.4570 

2.468 

0.44 

2.457 

2.447 

0.16 

2 

7 

1 

2.1354 

2.148 

.59 

1.866 

1.875 

.48 

3 

6 

2 

1.9826 

1.975 

.38 

1.674 

1.679 

.30 

4 

5 

3 

1.8577 

1.871 

.71 

1.578 

1.584 

.38 

5 

4 

4 

1.7861 

1.775 

.62 

1.744 

1.733 

.63 

6 

3 

5 

2.3877 

2.400 

.52 




7 

8 

2 

2.0395 

2.033 

.32 

1.775 

1.759 

.90 

8 

7 

3 

1.9295 

1.933 

.18 

1.689 

1.682 

.41 

9 

6 

4 

1.8367 

1.827 

.53 

1.676 

1.670 

.36 

10 

5 

5 

1.7861 

1.797 

.61 

1.862 

1.865 

.16 

11 

4 

6 

2.0147 

2.012 

.13 




Average 





0.46 



0.42 




'TABLE II - PRESSURE COEFFICIEHT POR EACH WING REGION 


typo of 
region 

Reference figure 

Pressure coefficient 

Special features 

A 


C P* S-[ H+P 2 + P 2.^ 

Affected only by supersonic 
leading edge. 

• 

B 


°P B §[ H + P 1 + p 2 ] 

Affected by subsonic and 
supersonic leading 
edges. 

C 


C P- M H+3?2 ] 

Affected by subsonic trail- 
ing and supersonic lead- 
ing edges. Kutta- 
Joukowskl condition 
Imposed, F^ m 0. 

D 


Cp= ^[h + Pi + Pt,i] 

Affected by supersonic and 
both subsonic leading 
edges. found for 

strip ad ’in same manner 
as F-^ for strip He. 

For F 1(l use Rj - || 

and measure 8^ j clock- 
wise at point A. 

E 


C P = ^[-H +P l +F l,l] 

Affeoted by supersonic and 
both subsonic leading 
edges. j determined 

as for region D with 

Rj “|f- 
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TABLE II - PRESSURE COEFFICIENT FOR EACH WING REGION - Concluded 


Type of 
region 

Reference figure 

Pressure coefficient 

Special features 

F 

/ 

C p - §• [-H + F 1(l ] 

Affected by one subsonic 
leading and one sub- 
sonic trailing edge. 
Fl^j determined as 
previously with Rj * a|=. 

Kutta-Jowkowski con- 
dition imposed on sub- 
sonic trailing edge. 

Fj_ * 0. 

G 


°» ■ f M 

Affected by both subsonic 
trailing edges. Kutta- 
Joukowski condition 
imposed. 

P 1 “ p l,l = °- 

H 

m 

°» • 5H 

Affected by both subsonic 
trailing edges. Kutta- 
Joukowski condition 
imposed. 

P 1 “ * 1,1 = °* 

i 



Shaded region is affected 
by externally Interact- 
ing flow fields. 
Solution not handled 
by graphical method. 
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Region influenced by 

supersonic leading edge 
Region influenced by 
subsonic leading edge 
Region influenced by 
subsonic trailing edge 


Supersonic leading edge 


Subsonic 
leading 
edge 


c 


-Supersonic 

trailing 

edge 


Figure 1. - Illustration of wing regions. 
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M Vj/Ug) 



rotation of terms In equation (1). 
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Figure 5. - Schematic hyperbolic integration grid. 
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(a) Bangs of F e from 0 to 2.1. 

Figure 8. - Variation of faotor Fg with { ) /c for various values of ©2 

at K * V2. 
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(o) Range of Fg from 0 to 8.8. 

Figure 8. - Concluded. Variation of factor Fg with ( v w -Vq ) /e for variou* 
valuea of 8g at K = Vs. 
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Figure 9. - Parameters in region A (v < 0). 
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, - Determination of C p for specific wing boundary 
at M • 


Figure XO 
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Figure 12. - Effect of curvature on factor F g at M » /^ 2 , 
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Two-dimension*! re Rim* 
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Figure 14..- Lift distribution for U mtfz on olrculsr-tlp wing swspt bsolt «t SO°. 
S& on aubaonio leading edge » oo, on aubsonio trailing edge ■ 0# 
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Figure 15. - Wing illustrating regions considered in table II 







